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A NEW CLASS OF GENERALIZED POWER LINDLEY DISTRIBUTION 
WITH APPLICATION TO LIFETIME DATA 


BRODERICK O. OLUYEDE, TIANTIAN YANG AND BOIKANYO MAKUBATE 


ABSTRACT. In this paper, a new class of generalized distribution called the Kumaraswamy 
Power Lindley (KPL) distribution is proposed and studied. This class of distributions con- 
tains the Kumaraswamy Lindley (KL), exponentiated power Lindley (EPL), power Lindley 
(PL), generalized or exponentiated Lindley (GL), and Lindley (L) distributions as special 
cases. Series expansion of the density is obtained. Statistical properties of this class of 
distributions, including hazard function, reverse hazard function, monotonicity property, 
shapes, moments, reliability, quantile function, mean deviations, Bonferroni and Lorenz 
curves, entropy and Fisher information are derived. Method of maximum likelihood is used 
to estimate the parameters of this new class of distributions. Finally, a real data example 
is discussed to illustrate the applicability of this class of distribution. 


1. INTRODUCTION 


Lindley [11] used a mixture of exponential and length-biased exponential distributions to 
illustrate the difference between fiducial and posterior distributions. This mixture is called 
the Lindley (L) distribution. Ghitany et al. [7| studied the statistical properties of the 
Lindley distribution. Sankaran [14] obtained and studied the Poisson-Lindley distribution. 
Jones [9] explored the background and genesis of the Kumaraswamy (Kum) distribution 
(Kumaraswamy |10]) and, more importantly, made clear some similarities and differences 
between the beta and Kum distributions. Among the advantages are: simple normalizing 
constant; the distribution and quantile functions have simple explicit formula which do not 
involve special functions; explicit formula for moments of order statistics and L-moments. 
However, compared to Kum distribution, the beta distribution has the following advantages: 
simpler formula for moments and moment generating function (mgf); a one-parameter sub- 
family of symmetric distributions; simpler moment estimation and more ways of generating 
the distribution via physical processes. Gupta and Kundu [8] provided a review and recent 
developments on the exponentiated exponential distribution. Cordeiro et al. [4] studied the 
Kumaraswamy Weibull (KW) distribution and applied it to failure time data. 

Motivated by the advantages of the generalized or exponentiated Lindley distribution with 
respect to having a hazard function that exhibits increasing, decreasing and bathtub shapes, 
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as well as the versatility and flexibility of the Kum distribution in modeling lifetime data, 
we propose and study a new class of distributions that inherit these very important and 
desirable properties, and also contains several sub-models with quite a number of shapes. 

In this article, we propose a new distribution, called Kumaraswamy Power Lindley (KPL) 
distribution. We discuss some structural properties of this distribution, derive the Fisher 
information matrix and estimate the parameters via the method of maximum likelihood. In 
section 2, we present some generalized Lindley distributions including the power Lindley dis- 
tribution, Kum distribution, Kum-G distribution, and the corresponding probability density 
functions (pdf). Section 3 contains results on the generalized and KPL distributions, includ- 
ing the hazard and reverse hazard functions, monotonicity property, and the sub-models. In 
section 4, we present the moment of the KPL distribution. Reliability and quantile func- 
tion are given in sections 5 and 6, respectively. Mean deviations are presented in section 
7. Section 8 contains results on Bonferroni and Lorenz curves. Measures of uncertainty, 
Fisher information and distribution of order statistics are presented in section 9. Maximum 
likelihood estimates of the model parameters and asymptotic confidence intervals are given 
in section 10. Section 11 contains an application of the proposed model to real data, followed 
by concluding remarks in section 12. 


2. SOME GENERALIZED LINDLEY DISTRIBUTIONS 


In this section, some recent generalizations of the Lindley distribution are given. The 
Kumaraswamy-G distribution is also presented. A useful series representation is given below. 
For |w| < 1 and b > 0 a real non-integer, we have the series representation 


(=u) t= cy (’ F e. 
j=0 
The one parameter cdf of the Lindley distribution [11] is given by 
LA ACRAS. 
dC ee 
The corresponding Lindley pdf is given by 
2 
e) m3) = AEn, 
Lindley distribution is a mixture of exponential and gamma distributions, that is f(x; A) = 
(1 — p) falx; A) + pfg(z; A) with p = D where falz; à) = GAM(2,A), and fg(z;4) = 
EXP(A). Now, let Yı and Y? be two independently gamma distributed random variables 


(1) Gr(z;à)=1 for x > 0, and A > 0. 


for z > 0, and A > 0. 


with parameters (o, A) and (a + 1, A), respectively. For y > 0, let X = Yi with probability 


es and X = Y; with probability 7—, then the pdf of X (see Zakerzadeh and Dolati [16]) 


is given b 
: i (Ar)? (a + yaje??? 
(A+ y)[(at+1) ' 





far; Q, à, y) = 
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for x > 0, à > 0, a > 0, y > 0. Note that when a = y = 1, we obtain the Lindley pdf given 
in equation (2). When y = 0 we have the gamma pdf with parameters a and A. If a = 1 
and y = 0 the resulting pdf is the exponential pdf with parameter A. 


2.1. Generalized Lindley Distributions. Nadarajah et al. [12] studied the mathematical 
and statistical properties of the generalized Lindley (GL) distribution. The cumulative 
distribution function (cdf) of the GL distribution is given by 


1+A+ Ax 
LA 


and the corresponding GL probability density function (pdf) is given by 


(3) Gaio X) = [1 - ew. 


qA? O 1+A+Àz 


(4) gaLlz; a, À) = x d +2) h LEA ea) | exp(—Az), 





for x > 0, A > 0, a > 0. This distribution is essentially the exponentiated Lindley dis- 
tribution. Zakerzadeh and Dolati [16] presented and studied another generalization of the 
Lindley distribution. These generalizations of the Lindley distribution are considered to be 
useful life distributions and are suitable for modeling data with different types of hazard rate 
functions: increasing, decreasing, bathtub and unimodal. These models constitute flexible 
family of distributions in terms of the varieties of shapes and hazard functions. 


2.2. Power Lindley Distribution. Ghitany et al. [6] presented results on a two-parameter 
Lindley distribution and referred to model as the power Lindley distribution. Considering 
the power transformation X = TV“, the cdf and pdf of the power Lindley (PL) distribution 
are given by 


ld4A-Az? ye 





(5) G pr (x; Q, A) = l= Tux € , 
and 

aa” Ka AA Wa ai WI Wa 
(6) gpL(£;a,A) NU 1 d n tc )v € , 


for x > 0, a > 0, A > 0, respectively. 


2.3. Kumaraswamy-Generalized Distribution. Kumaraswamy [10] in his paper pro- 
posed a two-parameter distribution (Kum distribution) defined in (0, 1). Its cdf and pdf are 
given by: 


F(r;a,b =1—(1—27)’, and f(rz;a,b) = abx^-! (1 — x2)", 


respectively, for x € (0,1), a > 0, b > 0. The parameters a and b are the shape parameters. 
Let G(x), be an arbitrary baseline cdf in the interval (0, 1). The Kum-G cdf F(r;a,b) and 
pdf f(z;a,b) are defined by 


(7) F(z;a,b) 21- (1 - [G(X)", 
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(8) f(z;a,b) = abg(zx)|G(x)]^ !(1— [G(zx)|?)" ! for a > 0,b > 0, 


dG 


respectively, where g(x) — dG) is the pdf corresponding to the baseline cdf G (x). 


3. KUMARASWAMY POWER LINDLEY DISTRIBUTION 


Now, with the choice of G(x) in the Kumaraswamy generalized distribution as the power 
Lindley distribution, we obtain the Kumaraswamy Power Lindley (KPL) distribution. The 
four-parameter KPL cdf and pdf are given by 


a ay b 
(9) Fkpr(z;a,à,a,b)=1— fı — h — DOE eec) | ; 
and 
(10) fkrr(z;a,A,a,b) = ab[Gpr(x) 1 — [Gez(x)]*]" ger(x) 
abad? 





= yb + se Ja?! exp( - Ax?) 


ee s 
— qb. ex^) 


14 
B: b 
fi- h- DM sw]. 


LA 


x 


(11) 


for x > 0, œa > 0, à > 0, a > 0, b > 0, respectively. Figure 1 illustrates some possible shapes 
of the pdf of the KPL distribution. 


3.1. Expansion of Density. In this section, the series expansion of the KPL pdf is pre- 
sented. When b > 0 is real non-integer, we use the series representation 
ajb—1 = pex] 
[1 — [Ger (x) = 5 (71) ; Jg» Js 


i=0 


where Gpr(r; A) = 1— 444“e™" | Tf a is an integer, from the above expansion and 





equation (10), we can write the KPL density as 


(12) fkpr(z;a,A,a,b) = abgpr(x pet uf 7 ) [G pr (z)] 0971 
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FIGURE 1. Plots of the pdf of KPL distribution for selected values of the parameters 








= ab ox (1+ 2°%)2°"! exp( -C Ax?) 
x yy (’ pi ) Í — CIE exp(—Ae*) us 
= ox (1+ z*)z*7! exp( Az?) 

(13) x > di h — ae exp(—Ae*) R : 


where the coefficients d; are 


YA 


d; = d;(a,b) = cya 2 ) 


and 9 "od, 1, for xr >0,a SOAS 0.0 0, b 0. 


If a is real non-integer, we can expand [G p; (a)]|*^?-! as follows 


[Gpz (z)]| 091 = ü-[- Geistes 
= Vy (" m D i- Gx (2)), 


j-0 d 


with 
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so that 
(14) (Gpz(a) eC = B Y cas Er zs = 3 (7) [C pr. 
j=0 r=0 
From equations (12) and (14), the KPL density can be rearranged in the form 
(15) fkprr(z;o,A,a,b) = gpr(x) >” aiea 
" i,j=0 r=0 
Ss (1 + 2*)2° 7 exp(- Ax?) 2 ` di jr 


,j20 r=0 


(16) ex^) i : 


where the coefficient d; jr is 


di, = di, (a, b) = (—1) tab i e: 5 ) e ) B 


and SU Bom dijr = 1, for x > 0, a > 0, A 7909 50, 5-0. Hence, Tor any real 
non-integer a, the KPL density is given by three (two infinite and one finite) weighted power 
series sums of the baseline cdf Gp;(z). By changing 5/7, Ss 2 Qm > s, in equation 
(16), we obtain 


oo 


fpx(a;o,d,a,b) = gpr(x) >. eilGpz (x) 


1,r—0 





= L+A+ Ax® ý 
17 x i |1 — ——— —Az? ; 
a? D a f= AE eec) 
where the coefficient c; is 


c; = ci(a,b) = cara z Herat Fi) —1), 


with 
WAA YA Aa X pr p E = ') H 


for x > 0, a > 0, A> 0, a > 0, b > 0, respectively. Note that the KPL density is given 
by three infinite weighted power series sums of the baseline distribution function G p; (xz). 
When b > 0 is an integer, the index 2 in the previous series representation stops at b — 1. 
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3.2. Some Sub-models of KPL Distribution. In this section, we present the sub-models 
of KPL distribution for selected values of the parameters a, a, and b. 
(1) b —1 


If b = 1, this is the exponentiated power Lindley (EPL) distribution with cdf and 
pdf given by 





Fgpr(z; œ, à, a) = [ — L eo^) ; 
and 
aan Mon » 
fepr(z; a, A, a) px epus Se) 
1-4 A Ax Ee 
i- TEX (Ar) 
for x >0,a>0, A> 0, a> 0, respectively. 
(2) a= 1 


If a — 1, the KPL cdf and pdf reduce to: 


1 a a 
TA TE E [HAPE opa) | 


and 


pone, 
14A 


1+A+ Az? xus 
E exp(—Ax | , 


for x > 0, a > 0, A> 0, b> 0, respectively. 
(3)a=b=1 


If a = b = 1, this is the power Lindley (PL) distribution given by equation (7). 
(4)a=1 


If a = 1, this is the Kum Lindley (KL) distribution with cdf and pdf given by 


1 ay b 
Fxr(z;A,a,b) 21— fı — [ — CC eC) | ; 


fkpr (x; Q, A, b) = Lee Ve T “exp(—Azx*) 


and 
aba? 
1+ 


LEA FAT "o 
- iat ae ex) 


LAHA AA 
x fı — h — CCS eC) | ; 


for x > 0, à > 0, a > 0, b > 0, respectively. 





fex(a; à, a, b) xd +2) exp(—Az) 
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(5)a=a=1 
If a =a = 1, then the KL cdf and pdf are given by 


LA FAT 


b 


Frola; A, 6) 21— | 


and 


bA? 
frels; à, b) = T (1+ a exo) | 


for x > 0, A > 0, b > 0, respectively. 

(6)a=b=1 
If a = b = 1, this is the generalized Lindley (GL) or exponentiated Lindley (EL) 
distribution (Nadarajah et al. [12]) given by equation (1). 

(D od bz 
Ifa = a =b = 1, this is the Lindley (L) distribution (Lindley [11]). The Lindley cdf 
and pdf are given by 


LA FAT 


b—1 
~ eC) | 





LA FAT 
Fa; A) == Tay CPC) 
and 
M? 
filz; à) = ix n + x)exp(—Azr), 


for x > 0, A > 0, respectively. 
(83)a=b=la=2 
If a=b=1,a = 2, then the KPL cdf and pdf reduce to: 
LAHA _\ 2 
1 c 





Fpy (zm; A) =]— 


and 
DAT 


1+. 





(1+ az?)me ^, 


fer(z; à) = 
for z > 0, A > 0, respectively. 
(9)a=la=2 
If a= 1,a = 2, then the KPL cdf and pdf reduces to: 





L+A+ Az? : 
PesrUmpASD)—e1- [ANE exe) , 
and 
2b)? 
frel; à, b) = 12: n + z?)z exp( C Az?) 
1+rA+ Az? oe 


for x > 0, A > 0, b > 0, respectively. 
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(10) b— 1,0 —2 
If b= 1,a = 2, then the EPL cdf and pdf are given by 





1 2 1 
fou h 7 A eA?) | 
and 
2a)? j $ 
fepL(x;à,a) = 1: 3 M PR )x exp( -Az^) 
1+A+Az? se 


for x > 0, A > 0, a > 0, respectively. 


3.3. Hazard and Reverse Hazard Functions. In this section, the hazard and reverse 
hazard functions of the KPL distribution are presented. Graphs of these functions for selected 
values of the parameters a, A, a, and b are also presented. The hazard and reverse hazard 
functions of the KPL distribution are given by 

fkpu la; Q, À, a, b) 

Frpz(x; a, 4, a, b) 





hu (aye; aso) = 





1 a a-—1 
,_itAtast ex^) 


a ay —1 
x fı- h- LT epo) | ; 


and 
freez(2; Q, À, a, b) 
Fi prm; Q, À, a, b) 


2 
= PAA (eat yit! exp(—Ar°) 





T o à, a,b) = 





1+ 
AMAT gd 
5 i- caer Wa a ] 
14A Az* Paine eae 
Fees eel 


for x >0,a>0,A>0,a>0, b > 0, respectively. The graphs of hazard function of KPL 
distribution are shown in Figure 2. These graphs show the variety of shapes for the KPL 
hazard function including bathtub, upside down bathtub, deceasing, and increasing hazard 
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rate functions. This attractive flexibility makes the KPL hazard rate function useful and de- 
sirable for non-monotonic empirical hazard behaviors that are more likely to be encountered 


in real life or practice. 
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FIGURE 2. Plots of the hazard function of KPL distribution for selected values 


of the parameters 


3.4. Monotonicity Property. In this section, we discuss the monotonicity properties of 
the KPL distribution. Let 


1 a 
(18) W(z) = Gp (5a, X) 21 — HEATA exp( Aa), 
14A 
then from equation (11), we can rewrite KPL pdf as 
bad? 
Tepe, à,a,b) = Tox + z*)x9-! exp(— Ax?) 


x [Wa -IWG)Y", 


for x > 0, a > 0, à > 0, a > 0, b > 0. It follows that 





2 
log fxpr(x) = log (=) + log (1+ 2?) + (a — 1) log x — Az? 


(19) + (a—1)log W(z) + (b — 1) log [1 — [W (a)]*], 
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and 
dlog fxpr(z) _ az" | «-1 E n" 
— Hec CO MERC UR AQT 
W'(x) [W (2)]W" (a) 
+ (a— Iw b) 1- Wr 
BB UL D zm ar! a-l 
i 1 Fao o g 


where W'(x) = SUA) = a(l + 2%) 2°! exp(—Ax®). 


Analysis: We know that x > 0, a > 0, à > 0, a > 0, and b > 0, so that 








; dW(z) Aa = 
= = 1 | T aes x " ? ` 
W'(x) d i3 z^)z^ "exp(-Ax^)» 0, Vx >0 
If x — 0, 
lobe à 
W(x) ET exp(—Az*) — 0. 
If z — oo, 
lc A-cAx 
Ea = 1 
W(z)-1 EF exp(—Az?) > 1, 
since 
14- A c Ax? raze! 


= 0. 





li e 
Exod (1 + A) exp(Ax?) 2d (1 + A) exp(Ax?)(Aoxe-1) 


Thus, W(x) is monotonically increasing from 0 to 1. Now, since 0 < W(x) < 1,0 < 


[W(z)? < 1,Va > 0, [W(z)! »0,Va»0,0«1-[W(z)]? < 1,Va » 0, and W'(z) > 0. 


Then we have T 0, and MEDI > 0. Also, from z > 0, we have x^^! > 0 and 


xri >O. 
IfA>l,a<1,a<1,andb > 1, we get SUR IDE < 0, since (1—A)axz* ^! — Aaz?^^1 < 0. 
In this case, fk pr(x;o, A, a,b) is monotonically decreasing for all x. 





If A € 1, fepr(2;0,A,a, b) could attain a maximum, a minimum or a point of inflection 


according to whether 


d? log fpr(x) 
dz? 


d? log fk pi (2) 
dz? 


d? log fgrrlz) 
dx? 











<0, >0,or =0, 


respectively. 
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3.5. Shape of Hazard Function. Note that if £ > oo, then HAPIT exp(—Az®) — 0. 





Also, 
1+A+ Ax® x a—1 oo E, 
bec ane cu ) = x ) 
1+A+ rA2x° zak 
g s prd Y ) 
LAs Ae E 

Tea = 1) IFA exp(—Ax^). 

Consequently, 
a9 ba AP 
= PEDEM (1+)! CDA exp ( Abg?) 
If x > 0, then 
2 

(22) fkprr(z;o,A,a,b) ~ abo (a4-1)a—1 


The cdf of KPL distribution is 


1+ a ay b 
Fkpr(z;o,A,a,b) 21— fı — [ — CHE aa | ; 





frx 590,0 5 0,.A 00.0 2907 b> 0: 


If x — oo, then 44 exp(—Axz?) > 0. Also, 


LA FAT ü > Ya 1+A+4+ Ax? ; 
T — a = — ————————————— — o 
LA alas ) >, @ LA op Ay ] 

LA Hare 


T= E N E a 
UA exp(—Azx°), 





so that 


1+A+ Ax* 


b 
Fkpr(z;o,A,a,0b) & 1— E poe exp(—\e*) 


Also, 


1T+tAt+ Ax? 
1 — Frpzr(a;a,A,a,b) S as 


b 
EF exp(—Ax ] 


b 
= a3 FA + Az*)' exp( —Abz?) 





UP NS 
(23) ~ Go R: 


If x — 0, then 2 exp(—Ax®) — 1, and 1 — HIT exp(—Az?) — 0, so that 





a 


1 a 
pO exp(—Az)| — 0. 


1+A 
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Thus, 


2 


1 pi d 
Fkpr (v; a, A, a,b) 1— fı —b h — L eC) | 


1+. 


= % h = CHO emot) 





1+. 


L1+A+ Ax d 
Rd 1 — —————(1- Ax? 
i| 1+. ( e) 


Aa? ay 2-2 a 
Ei pe ee 
[y 
b20 
24 ~w r”. 
(24) (1+ Ay 
The hazard functions of KPL is given by 


Jk pili; a, A, a, b) Jk pua; a, A, a, b) 
25 h ; , À, 1b — r^ = ? 
29) EEEE AU Fkpr(zy;o,A,a,b) | 1— Fkpr(z;o, A, a,b) 
for x > 0, a > 0, à > 0, a > 0, b > 0. If x > oo, with equations (21) and (23) in equation 


(25), we get 

















a bo Ae Lg HD! exp Abg yd + A)? 
a? >> exp(—Abax®) /(1 + A)? 


= bart. 





hy pr (2; a, A, a, b) $ 


If x — 0, with equations (22) and (24) in equation (25), we get 

aba geet (1 Jd aja 
1 — bM gae / (1 A)e 

aba?" gy (Dol 

(1 A)? 





hgpr (x; a, A, a, b) 





: bA2 aa 
since qi 27^ > 0, as z — 0. 


4. MOMENTS OF KPL DISTRIBUTION 


In this section, moments of the KPL distribution are presented. The following lemma is 
proved by using the result given by Nadarajah et al. [12]. 
Lemma 1 
Let 


xn) U exp(—qz) dz. 





kman = EX (mo) (9) (1) ave we. 
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2. If m is an integer, we have 





l k+l 
7 m — AN (UN (k+1\ (—1)'n*T (p +w +1) 
Kime) = 33 3. II 
Proof. (1) If m is non-integer, then 


h i Lin M epos] " > " A b (21) 


and 
K(m,n,p,q) = gem 


x [ z?(1-- z)(1-- n 4 nz) exp[- (nl + q)a] dz. 


Furthermore, l is an integer, so that 


(1 n-nz) = Y G (n -- nz)* = Y (0 +r)", 


k=0 


and 





x z?(1 -- z)**! exp[—(nl + q)a] dz. 


Now, k is an integer, so that 








and d 
Ud NEM BEE Ge ("21 
Y | sett 1 + q)a] de 
SUD Cel ee 5 8 


(26) 
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(2) If m is an integer, the index | in equation (26) stops at m — 1, so that 


kno) - YS ("21 (I) en) Con en 


1=0 k=0 w=0 

















The s^ moment of the KPL distribution, say ju, is given by 
TA = | x° freer le: a, A, a, b) dz. 
0 


Let b > 0 and a > 0 be real non-integer, then from equation (16), we obtain 


ee um x Se cl are ee ( 








i,j=0 r=0 
oo 1 a T 
x i (1+ z?)a** 7! exp(— Ax?) [ — CHE" emet) da 
0 1+ A 
abr? «E. d ey ee Dd sev 
ssj ijr 
punc que d 
(27) x ji (1+ z?^)z? exp( — Az?) h — CHO" aeo) dz“. 
0 LA 


Let y = x, then equation (27) changes to 


Fe wx >. Spite Li - ) d E ) G 


i,j=0 r=0 





e & LAHA d 
x / (1+ y)y* exp(—Ay) h — At am dy. 
0 


Applying Lemma 1, with m =r +1, n = À, p= %,q=4, we have 


ds DA » Spite s Ei = ) (’ a ) (7) 


0 
(28) x K(r+1,=,a). 











As r is integer, then r + 1 is integer, then the s‘” moment of the KPL is given by 


e BEEEC ICOMOS 


1$,7—0 r=0 1-0 k=0 w=0 
(-1) HT (2 +w + 1) 
(1 Alia tw=k-1(] plijeterr 


Note here that we have considered the case when b > 0 and a > 0 are non-integer, however 











(29) 





the other cases can be similarly derived. 
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5. RELIABILITY 


In reliability and related areas, the stress-strength model describes the life of a component 
with random strength X, that is subjected to a random stress Y. The component will fail at 
the instant that the applied stress exceeds the strength, and the component will function sat- 
isfactorily whenever X > Y. We derive R = P(X > Y), a measure of component reliability, 
when X and Y have independent K PL(a1, 1, 41,61) and K PL (ag, A2, a2, b2) distributions, 
respectively. Note from equations (9) and (11) that 





























R = P(X>Y) 
= | fx(x; a1, A1, a1, b1) Fy (2; a2, A2, a2, b2) dx 
0 
id a1b104A7 ESSI 
= ———=(1+a4™%)r2™~ exp(—A,2™ 
[PER es) essa) 
1 + Al + Age M a1—1 
h — ES exp(-Aix ] 
1 ai ay» bı—1 
pe pitt 1^ 
= aj1b10, A? Ly 
—— —— — (4 + z?3 )z?* ^^ exp(—A,2™ 
] PERO ate eoa) 
14 AG + azi m a1—1 
[ - ESS exp(-Aix | 
1 ai ay» 51—1 
mae) 
as à2 b2 
(30) x fı — [ E = = exp(—dxe") | dz. 





Applying the series expansions 


1+ M1 + MEO HW E ay, — 1 k 
i ai = 
a ee SE 
(— 


1) Am 4": exp( — A; kx?!) 
(Pe Aa) 








x 





(eS eee) ex CT 


l,p=0 n=0 
(=1) ? Ma"? exp(-Aipa?) 
(14- A1)" 
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16 + Aage? 2: à2 x b2 
fı — h = D exp(—Aax | j 


e = SO) 


g,t=0 h= 








and substituting equations (31), (32), and (33) into equation (30), we get 


pda x GT LG 



































x (sq) ttp ymn? fo + eMart al ex (x (k E La) dz 
(1+ Artt Ju P 1 p 
oo k p t 
^ —^v^wv[1-— 1 k bi = 1 al p bo 
um Ce CUR) 
dE k /NmJN I n AN, 
aoq\ (t (1) ar tnt Nk 
t h) (14 A,)m*tnH (1 + A9)^ 
(34) x T (1 + gory iret aithaz-1 oxn(— ACE p Daa z Asta?) dex 
0 
Note that 
| (1 + £21) amini Da oxn(— A (Kc +pt 1)x?!) dz 
0 
1 oo 
= — glmtnjar exp(—Aj(k +p+ Dar) dze! 
O1 Jo 
1 oo 
ox gp lmtnt or exp(—Ai(k +p+ Dai) dz 
Qı Jo 
1 oo 
= — | y™™exp(—Ai(k +p + 1)y) dy 
O1 Jo 
1 T m+n+1 
up exp(—Ai(k + p + 1)y) dy 
1 Jo 
= T(mtn+])  , — T(mtn-2) 
oa [X (1 k pier ! oa [A1 (1 + k + p) ++? 
(35) = (m+n)! p ment 
oipa rk p) | X kp" 
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oo (—-1)hftiaie2 


where we let y = x^'. Also, by using exp(—Astz?^?) = $775 E , we get 


(1 + go1)g mt rDoihoz-1 exp(— i(k + p+ 1)z?* — Agta?) dz 


























i 

= f glmint ortho! oxp(—\1(k + p 4- 1l)a™ — Agte™) dx 
0 
i 









































+ gm Da tha exn(— i(k + pt 1)z?! — Agta™) dz 
= —1 Piu SE | fi fi L4 
e y» ( Post | gimtntlort(h+i)az ! exp(—A4(Kk p 1)z?1) dz 
l: 0 
i—0 
= —1 DATA » | fi fi L4 
2» , | ge a Lb ep Ea d 
i 0 
i—0 
c (-DA f 
= > EA z 2 P UE ERO Ga lita) ee exp(—Ai(Kk + p+ 1)xz?!) da™ 
io 2:01 0 
= (-1)XBé f° 
| ` ( » 2 | glmtnt tor Hijaz exp(—A1(k +pt Da?) dz? 
i—0 2:0 0 
iU 2:01 0 
| y ( 3 2 f y +n+1+(h+i) 2 exp( —Ai(k + p + 1)y) dy 
pum 2:01 0 
_ 3 (—1)'d3t! 
0 ila, 
(36) ZR 
[A (ine | acne [A (14 k4 jperexer um ! 


where we let y = x. Substituting equations (35) and (36) into equation (34) to get 


me PIG) 


(—1)F**»A (m+n)! | , mtn+l | 
(1+ A4)netn1(1 + k + p)mteH | Ai(1 9 Kk 4 p) 
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0101 k m l p n q 
k,l,p,q,t,i=0 m=0 n=0 h=0 


—(h+i) 22 ns 
1 ang 


ird (1) (—1) bte tH y 
t h (1+ A,)m (1 + Ag)^(1 + k H p) iH 092g 
D mspdsp2 s (eat) 


1+k+p 














6. QUANTILE FUNCTION 


The quantile function, say Q(p), is defined by F(Q(p)) = p. Now, from the cdf of the KPL 
distribution, we have 


b 


Esc MOPE epale] , - 


LA 





Fxpx(Q(p)) = 1- fı - h ii 


and we can obtain Q(p) as the root of the following equation 


i 
a 


IA HALO (pe exp(—A[Q(p)]*) = E a »* Eu 


(97) 1+A 





for 0 < p < 1. Substituting Z(p) = —(1 + A + A[Q(p)]^), we can rewrite equation (37) as 


ZD exp(1 ++ 20) = [i-a - 9]* -a 


so that 
zoezo) = (1+ Xe C1 - x { [1 = 0 - 9 te it. 


for 0 < p < 1. As the defining equation for Lambert W function W (x) is x = W (x) exp(W(x)), 
we get 


Zip) = W (a +y- { [1 - a -p)]*- p. 


for 0 < p < 1. Then, we obtain 





Q(p) — 3 ; 


fr0<p< 1. 
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7. MEAN DEVIATIONS 


The mean deviation about the mean and the mean deviation about the median are defined 
by 


aX) = f le- mfa) de, and BX) = f le- MIfi)ds 


respectively, where u = E(X), and M = Median(X) denotes the median. The measures 
ôı(X) and ô2(X) can be calculated as follows: 


(38) AQ) = Fu) -2042 | eo 
and 
(39) BX) = -u+2 | afir)ds 


respectively. By using the moments for KPL distribution and the results in Lemma 2 
(Nadarajah et al. [12]), we can calculate equations (38) and (39). Note that 


l+tnt+naz 


m-—1 
Een amna) exp(—qz) dz, 


K(m,n,p,q) = [va + 2) h- 


and 


ltnt+nez 


m-—1 
lia apona) exp(—qz) dz. 


Limnpat) = f eats) fi- 
t 


We consider the case when b is real non-integer and a is non-integer. From equation (28) 
and Lemma 2 (Nadarajah et al. [12]), we know that 


EEE NIG 








H = 
1 1,9=0 r=0 
1 
(40) x K PE LA za). 
a 
i ab? C cotes üt Lege tbe) 
«às = PY Y cose( 197g: 
J TrA a J t E 
(41) x L (+ -1,A, D : 
Q 


and 





M 


[aene = Eiet- ON) 
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so that 


5,(X) 2uF(u) - 2n + 1 a =o) c pu $ li ) (’ i ) G) 


Panda’ 
a 





x 


and 





Bera a x Spit LI - ) ( : n () 
x L(r+1, à: > AM). 


Note here that we have considered the case when a and b are non-integer, however the 


other cases can be similarly derived. 


8. BONFERRONI AND LORENZ CURVES 


Bonferroni and Lorenz curves are defined by 
q 1 f 
Boe) = — f xf(x)dz, and up - f ap ay de 
respectively, where u = E(X), and q = F-!(p). Now, we obtain Bonferroni and Lorenz 
curves for KPL distribution as follows: If b is real non-integer and a is non-integer, then 


from equations (40) and (41), we have 


Peroa = S Eae (này zn) 





Bor = rT ad 


pater IIE 
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and 
ab US yer a(1+i)—1\ /b — 1X fj 
rediere X Qe EU 
i,j=0 r=0 
x don Eaa), 
respectively. 


Note here that we have considered the case when a and b are non-integer, however the 
other cases can be similarly derived. 


9. ORDER STATISTICS, MEASURES OF UNCERTAINTY, AND INFORMATION 


In this section, the distribution of the kt? order statistic, measures of uncertainty, and 
information for the KPL distribution are presented. The concept of entropy plays a vital role 
in information theory. The entropy of a random variable is defined in terms of its probability 
distribution and can be shown to be a good measure of randomness or uncertainty. 


9.1. Distribution of Order Statistics. Suppose that X,,---,X, is a random sample of 
size n from a continuous pdf, f(x). Let X44, < Xan < +++ < X44, denote the corresponding 
order statistics. If X,,--- , Xn is a random sample from KPL distribution, it follows from 
the equations (9) and (11) that the pdf of the kt? order statistic, say Y, = Xin, is given by 





hk) = an Wake ax "ye Y Fez (ye) Fee (Ye) 


n k—1-H j 


_ aboXnY + yg)yg  exp( Aug) APON Y >” 
(1+A)(k- 1)\(n — k)! l p qi, j=0 r 


ey oe a 


(esp PUR [W (yp) t", 





J 











x 


x 





where W (yk) = Gpr(Yk; a, A) 21— EE exp(—Ayf). The corresponding cdf of Y; is 


Hw) = $ 3 200) (3029 E owa. 


The st” 
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moment of the kt? order statistic Y; from KPL distribution is obtained as follows: 


If b is real non-integer and a is non-integer, then 


E(Y;) 


x 


x 


f EAEE 
0 


n—kk-—1-H ow 
aban! (— SI 


EEA BPS 


Clara t e) 


ni (1 -- yu exp( - Ay) W (yx) dg 
0 








p. 









































pee Edi ires 
Sc n—k aro Men 0 r=0 p 

bp\ (a(1--$) — 1V (b — 1X (3 

q J 1 T 

N SD 7 LEAH aye alg RS 
a (1+ Ye) Uk PO AY) Í = = oC) dy; 
abNn\(—1)tetatitste A >. on (" = J ( aries : 
UHA — 1)!(n — k)! — en Nerei es l p 

bp\ a1 +i) — 1X f/b-1\ fj 

q J 1 T 

m T Ice ee De KIT 
| (1 + t)te exp(—At) |1 — ESA eC) dt 
aber (Ie rate TER LIE oo mE (eil 
(1 A)(k — 1)l(n — k)! —— eM l X p 








ie te ee b— 

q 

ORT 
a 


where we let t = yẹ. 


Note here that we have considered the case when a and b are non-integer, however the 


other cases can be similarly derived. 


9.2. Renyi Entropy. Renyi entropy [13] is an extension of Shannon entropy. Renyi entropy 


is defined to be H,(fxpr(x)) = 


log( Je” fè pr (ia,r,a,b) dr) 
1—y 





H, (fpr (a3, À, a,b)) = , where y > 0, 
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and y Æ 1. Renyi entropy tends to Shannon entropy as y > 1. Now, 


e (aba? ? 
/ ficpy(@) dx zx (=) 


/ Gra Sepe 
0 


x [W(z)^- - [W e de. 





x 








(42) 
Note that 
weg = fi- aa 
Se us * o C) 4 29) 
(43) A oul JA ? 255 PA ve ) exp(—Akz?), 
and 
nwon = Yee 
(44) PE WA ue = exp(—Anz*). 


Substituting equations (43) and (44) into equation (42), we get 
Se abo? V! E ERE 
frpp(@)dz = ( ) ` » 
f 1 + À k,m,n=0 j=0 t=0 
ay — YN f KV (by — 7N (am) (mN (-1)Ft™ ett 
k j m no Xt (1+ Ajit” 


ri (1+ zo) ttg (079 exp( — A(4 + k + n)x?) da. 
0 








x 


(45) 


By using (1+ ze) t+ = Y 7*9 (7*7*5) zi in equation (45), we get 


/ (1 + rep Hrd exp(—A(y + k + n)z?) dx 
0 


























l1 HEN [S oq 

=% i J ) , TUKI Da exp( — Ay + k + n)z?) da? 
ao i r 
Lx IHA pee sy 

= — ü 7 ) f gia exp(—A(y + k + n)y) dy 
Cr i i 

map! (7 j i Dow 
am O G+ ken 
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where we let y = z^. Now, equation (45) simplifies to 
oo k n ay ^j k 
Ved E 
Eu D) 
k,m,n,i—0 j—0 t=0 
by —7\ fam\ (nV (y 4-3 t 
m n t i 
(-1) mi p(G +y+ 7) 


(46) x ee N 
(1 + AJ trt 71753 thra (y+ k-c-n)'rt 


ri fepi(*) dx 


x 











Consequently, Renyi entropy for KPL distribution reduces to : 





Hy(fxeu(e)) = z boglab) a) 


tref X: ESOO 





k,m,n,i—0 j= 


(—1) tmn (i +y + ty) 
Gyan i Sy 4 k p nits J’ 








for y > 0, and y Z 1. 
Note here that we have considered the case when a and b are non-integer, however the 
other cases can be similarly derived. 


9.3. s-Entropy. The s-entropy for KPL distribution is defined by 


a = Jo. Fess (tte Ava, ba 3E $5 1. 8 9.0, 


Hifp; Q, A, a, b)) D Bs log / (X)] if es 


Consequently, if s Z 1,5 > 0, then from equation (46), we have 


Haa = —L.- 69 Y 303 es JC) 


k,m,n,i—0 j—0 t=0 


ec 
(—1)& m^ T (i 4. s + is) 


a 


(1+ )sthtn yisi- 2" (s +k+ njiti 





x 








If s = 1, then s-entropy is Shannon entropy. 
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9.4. Fisher Information Matrix. This section presents a measure for the amount of in- 
formation. This information measure can be used to obtain bounds on the variance of 
estimators, and as well as approximate the sampling distribution of an estimator obtained 
from a large sample. Furthermore, it can used to obtain an approximate confidence interval 
in case of large sample. 

Let X be a random variable (rv) with the KPL pdf fkpL(: ©), where © = (61, 02, 03,04)? = 
(a, 4, a, b)? . Then Fisher information matrix (FIM) is the 4 x 4 symmetric matrix with ele- 


ments: 





L;(9) = Eo | eet ©)) Olog(f pr (X; e» | 


00; 00; 
If the density fjpr(.; ©) has second derivative for all 1 and j, then an alternative expression 
for L;(0) is 





(47) L;(9) = -Eo É log(frrz(X; 2 


00,00; 

For the KPL distribution, all second derivatives exist, therefore the formula above is appro- 
priate and, most importantly significantly simplifies the computations. The elements of the 
observed information matrix of the KPL distribution are given in Appendix A. 


10. MAXIMUM LIKELIHOOD ESTIMATORS 


In this section, the maximum likelihood estimates (MLEs) of the parameters a, A, a, 
and b of the KPL distribution are presented. If 71,--- ,£n is a random sample from KPL 
distribution, then the log-likelihood function is given by 


2 TL TL 
log(L(o,A,a,0)) = nlog (=>) + X log(1 +r) spear) > log Ti 
i=l i=l 





E 2m t (a— 1) 2 log (W(z;)) 


(48) 3" pest) 2r log — [W (z;)]*], 





14AAz* " 
where W(z;) 21— ~ exp(—Az?). 
The partial derivatives of log L(a, A, a, b) with respect to the parameters a, A, a, and b 


are: 








= E+ Sloan) x a c) PCIE DOT. 


O log L(a, A, a,b 
Ob 





) = R4 Slogi - Weal 
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O log L(a, A, a,b) n T r7 log x; | z L o 


























E e Wla) 
^ [W (z)]* (OW (z;)/8o) 
+ a(1- b) x 1— Wy , 
and 
OlogL(a,A,a,b) |  m(24 2) 7 a ^. OW (7;)/0A 
OA = NEN) D Tp) 2. Win) 
24 IW (2,)] 871 OW (2) /8X 
+ a(1-— b) x [W a woe ) 
where 
OW (xi) _ A?(1 + z?)z? log x; exp( — Az?) 
da 1+A i 
and 
OW (a) 1+2A4+ Az 1 " Aa 
od | Jap ad d Ta EXPE, 


When all the parameters are unknown, numerical methods must be used to obtain es- 
timates of the model parameters since the system does not admit any explicit solution, 
therefore the MLE (4, A, à, b) of (a, A, a, b) can be obtained only by means of numerical pro- 
cedures. The MLEs of the parameters, denoted by © is obtained by solving the nonlinear e- 








: OlogL OlogL OloggL OlogL AT _ g : 
quation ( acu rac um ee )* = 0, using a numerical method such as Women arem 
procedure. The Fisher information matrix given by KO) = [Ig,,Jaxa = E(—2,98L) 
1005 


i,j = 1,2,3,4, can be numerically obtained by MATLAB or MAPLE software. The to- 
tal Fisher information matrix I (©) = nI(©) can be approximated by 
0? log L 
06,00; 





(49) J (Ô) « | , 1,j-1,2,3,4. 





IPI 4X4 
For real data, the matrix given in equation (49) is obtained after the convergence of the 
Newton-Raphson procedure in MATLAB or R software. 


10.1. Asymptotic Confidence Intervals. In this section, we present the asymptotic con- 
fidence intervals for the parameters of the KPL distribution. The expectations in the FIM 
can be obtained numerically. Let Ô = (à, À, à, b)T be the MLE of © = (o, A, a, b). Under 
the conditions that the parameters are in the interior of the parameter space, but not on the 
boundary, the asymptotic distribution of /n(8 — ©) is N4(0, I^! (9)). 

The multivariate normal distribution with mean vector (0, 0,0, 0)7 and covariance matrix 


I !(0) can be used to construct confidence intervals for the model parameters. That is, the 
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approximate 100(1 — 7)% two-sided confidence intervals for a, A, a and b are given by: 


â+ Za), A+ ZaVIX(O), GcZiyIL4(8) and b+ Zy Ip (8), 


nA A 


respectively, where I71(8), I1 (8), I71 (8), and I5; (8) are the diagonal elements of I7! (8) = 


^aa 


(nI(8))7!, and Z z is the upper ath percentile of a standard normal distribution. 


11. APPLICATION 


In this section, application of the KPL distribution including the estimation of the param- 
eters via the method of maximum likelihood and likelihood ratio (LR) test for comparison of 
the KPL distribution with its sub-models for given sets of data are presented. The examples 
illustrate the flexibility of the KPL distribution in contrast to other models including the 
Kumaraswamy Lindley (KL), power Lindley (PL), GL, L, Kumaraswamy Weibull (KW), 
Weibull (W), and gamma (GAM) distributions for data modeling. 

The MLEs of the KPL parameters a, A, a, and b are computed by maximizing the ob- 
jective function via the subroutine NLMIXED in SAS. The estimated values of the parame- 
ters (standard error in parenthesis), -2log-likelihood statistic, Akaike Information Criterion, 
AIC = 2p — 2log(L), Bayesian Information Criterion, BIC = plog(n) — 21og( L), and Con- 
sistent Akaike Information Criterion, AICC = AIC + geet where L = L(Ó) is the value 
of the likelihood function evaluated at the parameter estimates, n is the number of obser- 
vations, and p is the number of estimated parameters are presented in Table 3. The KPL 
distribution is fitted to the data sets and these fits are compared to the fits using the KL, 
PL, GL, L, W, and GAM distributions. 

We can use the LR test to compare the fit of the KPL distribution with its sub-models for 
a given data set. For example, to test a = b = 1, the LR statistic is w = 2[In(L(4, A, â, b)) — 
In(L(à, A, 1, 1))], where â, À, à, and 6, are the unrestricted estimates, and @, and À are the 
restricted estimates. The LR test rejects the null hypothesis if w > x where x denote the 
upper 100d% point of the y? distribution with 2 degrees of freedom. 

Specifically, we consider a maintainance data set. The set of data is the maintenance data 
with 46 observations reported on active repair times (hours) for an airborne communication 
transceiver discussed by Alven [1], Chhikara and Folks [3], and Dimitrakopoulou et al. [5]. 
The data set is given in Table 1. The MLEs of the parameters with standard errors in 
parenthesis and the values of the statistics (-2In(L), AIC, AICC and BIC) are given in Table 
3. The starting points of the iterative processes for the data sets for the K PL(o, A, a, b) 
distribution are (1, 0.115, 0.026, 0.1). 

Probability plots (Chambers et al [2]) consists of plots of the observed probabilities against 
the probabilities predicted by the fitted model are also presented in Figures 3 and 4. For 
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0.2 0.3 05 0.5 05 0.5 06 0.6 0.7 0.7 0.7 
0.8 08 10 10 LO 10 11 13 15 15 15 
15 20 20 22 25 2.7 30 30 3.3 33 4.0 
40 45 47 50 54 54 7.0 7.5 88 9.0 10.3 
220 245 - - - - - - - = E 
TABLE 1. Maintenance Data [1], [3], [5] 











Data set Model a À a b —2ln(L) AIC AICC BIC SS 
I (n=46) KPL(o,A,a,b) 0.6679 4.4222 | 5.6016 . 0.1488 200.6 208.6 209.6 215.9 0.06484707 











(0.05909) (0.8266) (2.7768) (0.02612) 
KL(1, A, a,b) 1 1.6091 1.0576 0.2029 2127 2187 2192 2242 04270972 
(0.03375) (0.3385) (0.03180) 
PL(o,A1,1) 0.7581 0.6757 1 1 2100 2140 214.3 2177 0.1185674 
(0.07424) (0.1016) 
GL(1, A,a,1) 1 3677 0.6643 1 215.7 219.7 220.0 223.4 0.2635925 
06442) (0.1352) 
L(1,4,1,1) 1 4664 1 220.0 222.0 222.1 223.8 0.5676042 
04990) 
Wa,d,1,1) 0.8086 0.2949 1 1 2089 2129 213.2 216.6 0.1156807 
(0.09576) (0.05138) 
GAM(a,A) 0.9323 02585 - : 2009 2139 2141 217.5 0.1716121 
(0.1701) (0.06150) 














TABLE 2. Parameters Estimates, Log-likelihood, AIC, AICC, BIC, and SS 


the KPL distribution, we plotted for example, 


at 





"E. 1+A+Aye m 
(50) Fly â, à, â,b)=1-— {1-— |1— STAT Ak pl ; 
LA 
against ee k = 1,2,---,n, where y, are the ordered values of the observed data. A 


measure of closeness of the plot to the diagonal line given by the sum of squares 


n 


2 
2255 LED b- E23] | 
was calculated for each plot. The plot with the smallest SS corresponds to the model with 
points that are closer to the diagonal line. The KPL model performs very well in this regard. 
For the maintenance data, the LR statistics for the test of the hypotheses Ho: K L(1, A, a, b) 
against Ha : K PL(a,A,a,b), Ho: PL(a, 4,1, 1) against Ha : KPL(o,A,a,b), Ho : GL(1, 4, a, 1) 
against Ha : KPL(a,,a,b), and Ho : L(1,4,1,1) against Ha : KPL(o,A,a,b) are 12.1 
(p — value = 5.04 x 1074 < 0.001), 9.4 (p — value = 9.095 x 10^? < 0.01), 15.1 (p — value = 
5.2611 x 1074 < 0.001), and 19.4 (p — value = 2.2597 x 107“ < 0.001), respectively. Con- 
sequently, we reject the null hypothesis in favor of the KPL distribution and conclude that 
the KPL distribution is significantly better than the KL, PL, GL, and L distributions based 


30 OLUYEDE, YANG AND MAKUBATE 


on the LR statistic. The KPL distribution is also better than the Weibull and Gamma dis- 
tributions based on the values of the statistics AIC, AICC and BIC. The plots of the fitted 
KPL distribution and sub-models are shown in Figure 3. 


Density Fit Probability Plot 








O 
m 
























Observed Probability 
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x Predicted Probability 


FIGURE 3. Fitted densities and probability plots of KPL distribution and 
sub-models for maintenance data 


Based on the values of these statistics, we conclude that the KPL distribution provides 
a better fit than the KL, PL, GL, L, and GAM distributions. For the maintenance data, 
KPL distribution is far better than its sub-models, and a pretty good competitor to the KW 
distribution. The KPL model can provide better fits than other common multi parameter 
lifetime models. 


12. CONCLUDING REMARKS 


A new class of generalized Lindley distribution referred to as Kumaraswamy power Lindley 
(KPL) distribution with flexible and desirable properties is proposed. Properties of the 
KPL distribution and sub-distributions were presented. The pdf, cdf, moments, hazard 
function, reverse hazard function, reliability, quantile function, mean deviations, Bonferroni 
and Lorenz curves were presented. Entropy measures including Renyi entropy, s- entropy 
as well as Fisher information matrix for KPL distribution were also derived. Estimate of 
the model parameters via the method of maximum likelihood obtained and application to 
illustrate the usefulness of the model to real data given. 
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APPENDIX APPENDIX A FIM For KPL DISTRIBUTION 


Let l = L(a, 4, a,b), and W(x) = Gpr(z;a, A) = 1— EE exp(—Az?). Elements of the 
observed information matrix of the KPL distribution are given by 
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Note here that we have considered the case when b > 0 and a > 0 are non-integer, however 
the other cases can be similarly derived. 
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